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Anharmonicity changes the solid solubility of an alloy at high temperatures 
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We have developed a method to accurately and efficiently determine the vibrational free energy 
as a function of temperature and volume for substitutional alloys from first principles. Taking 
Tii-^ARN alloy as a model system, we calculate the isostructural phase diagram by finding the 
global minimum of the free energy, corresponding to the true equilibrium state of the system. We 
demonstrate that the anharmonic contribution and temperature dependence of the mixing enthalpy 
have a decisive impact on the calculated phase diagram of a Tii_ x AGN alloy, lowering the maxi¬ 
mum temperature for the miscibility gap from 6560 K to 2860 K. Our local chemical composition 
measurements on thermally aged Tio.5Alo.5N alloys agree with the calculated phase diagram. 


When discussing the solubility of an alloy, the config¬ 
urational entropy is always taken into account, but the 
effect of temperature associated with lattice vibrations 
is often neglected. It has been experimentally shown [1- 
3] that the vibrational and configurational entropies are 
comparable in the cases of fee Ni 3 Al and Cu 3 Au, and bcc 
Fe 3 Al, and theoretical studies draw the same conclusion 
[4-8] — that lattice vibrations cannot be ignored. High 
quality phonon spectra can be computed in the frame¬ 
work of density functional perturbation theory (DFPT) 
[9] and the small displacement method [10], but the in¬ 
troduction of substitutional disorder as in an alloy causes 
the cost of such calculations to escalate rapidly. In this 
work, we propose a method of computing the vibrational 
free energy of a configurationally disordered solid based 
on the temperature dependent effective potential method 
(TDEP) [11, 12], which has an efficiency comparable to 
the state-of-the-art methods that only apply to ordered 
solids. Moreover, the method has the advantage of tak¬ 
ing into account anharmonicity of the lattice vibrations 
and therefore remains valid at temperatures for which 
the quasi-harmonic approximation breaks down. 

We demonstrate the accuracy of our technique in a 
study of decomposition thermodynamics of TR,AR_ X N 
alloys [13, 14], a system for which lattice vibrations un¬ 
derpin an unusual and technologically useful isostruc¬ 
tural decomposition [15]. Metastable TGAR_ X N coat¬ 
ings are ideal for use in the manufacture of cutting tools 
due to their characteristic age hardening during use. 
Metastable cubic TiAIN undergoes spinodal decomposi¬ 
tion to form nano-scale domains of cubic TiN and AIN, 
through which extra stress is required to propagate dislo¬ 
cations [13, 16, 17]. Remarkably, calculated values of the 
maximum temperature for the miscibility gap vary be¬ 
tween approximately 6050 and 9000 K [18, 19], depend¬ 


ing on the methodological details, and as low as 3790 K 
within the Debye-Griineisen approximation [20]. These 
are well above the dissociation temperatures of TiN and 
AIN, however, but cutting tools may reach temperatures 
of up to 1300 K [21], at which point vibrations could 
be of considerable importance and are subject to the ef¬ 
fects of thermal expansion and anharmonicity. We note 
that existing studies either neglect or use an incomplete 
description of the vibrational contribution to the free en¬ 
ergy, since the methodological challenge and computa¬ 
tional efforts required to calculate the phonon spectra 
of a substitutionally disordered solid using ab initio ap¬ 
proaches are considerable. 

In this Letter, we propose a computationally tractable 
method for the treatment of vibrational free energy of 
a random alloy, and use it to perform accurate first- 
principles calculations of the vibrational free energy of 
B1 TR_ X A1 X N alloy, our model system. The theoreti¬ 
cal miscibility gap has a maximum temperature of 2860 
K, and the solubility limit of A1 in TiN at intermediate 
temperatures is increased in comparison to calculations 
which neglect the effect of lattice vibrations. Our method 
employs the TDEP method to compute the vibrational 
contribution to the free energy. When constructing the 
phase diagram we minimize the Gibbs free energy [22] to 
obtain the stable alloy compositions at equilibrium. We 
perform atom probe tomography experiments, exploiting 
the high spatial and chemical resolutions to depict the de¬ 
composition of a supersaturated solid solution of TiAIN 
with an alloy composition inside the miscibility gap, thus 
verifying the predicted phase diagram, and demonstrat¬ 
ing the importance of lattice vibrations. 

Phase stability at constant temperature and volume is 
determined by the Helmholtz free energy F , which can 
be expressed as the sum of the free energy of a model 
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system F m and a correction term based on the Kirkwood 
coupling theorem [23], where U m and U are potential 
energies of model and real system respectively. 

F = F m + [ ( U-Um) x d \. (1) 

Jo 
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Eq. 1 is formally exact for phase equilibria in alloys at 
relevant temperatures. In practice, calculating A F can 
be very difficult [24] , so our strategy will be to choose the 
model system with the smallest possible A F. Ideally, it 
should be within the error bars of ab initio calculations. 

Our model system consists of atoms distributed on the 
sites of a special quasirandom structure [25] (SQS), al¬ 
though any other supercell approximation of the com¬ 
pletely random alloy could be used. These atoms, how¬ 
ever, interact with effective force constants which have 
full symmetry of the underlying crystal lattice. Lattice 
vibrations of the model system are described by tempera¬ 
ture dependent effective potential (TDEP) model Hamil¬ 
tonian [11, 12]: 

^ = + (2) 
i 1 ij 

where p; and are the momentum and displacement of 

= eff 

atom i, are the second order effective force constant 
tensors, which relate the displacement of atom j to the 
force f) exerted on atom i. The TDEP method allows us 
to fit parameters of the model Hamiltonian to results of 
ab initio molecular dynamics simulations carried out for 
the real system of interest. 

In the model Hamiltonian (Eq. 2), the translational 
and spatial symmetry of the underlying crystal lattice 
are imposed by treating the components of the alloy as 
symmetry equivalent, and the real interactions present in 
DFT are mapped to effective interactions expressed by 

= eff 

the interatomic effective alloy force constants <f > , J corre¬ 
sponding to the underlying crystal lattice, B1 in case of 
Tii_ x Al x N. Imposing the symmetry reduces the number 

= eff 

of independent components in ^ ■ considerably, mak¬ 
ing the method numerically efficient [12] thus we call our 
approach the symmetry imposed force constants TDEP 
method (SIFC-TDEP). In our approach, the forces are 
obtained from calculations of the real disordered system 
as described by the SQS and only then mapped onto 
effective force constants in a system with randomly dis¬ 
tributed atomic species having their respective correct 
atomic masses. Note that the accuracy of the free en¬ 
ergy calculations can always be improved, for example, 
by using a more complex model Hamiltonian. This can 
be achieved using an expansion of the pairwise force con¬ 
stants AB, BC, AC... resulting in different pairs for dif¬ 
ferent pairs of species, although this is not necessary for 
our model system. 



FIG. 1. (Color online) Comparison of harmonic and anhar- 
monic vibrational Gibbs free energies as a function of temper¬ 
ature. Harmonic free energies were calculated via the TDEP 
(SIFC-TDEP-HA) and small displacement (SMD-HA) meth¬ 
ods. The free energy including anharmonic effects via the 
terms AUo and g(oj,T) (SIFC-TDEP) is compared with a 
single thermodynamic integration data point. 


In order to find the effective force constant matrix that 
best represents the Born-Oppenheimer potential energy 
surface, we minimize the difference in forces between the 
model system and the SQS model of a real alloy, comput¬ 
ing the latter by means of ab initio molecular dynamics 
(see supplementary material). 

We calculate the internal energy via Density functional 
theory (DFT), with the electronic contribution to the 
entropy given by the Merrnin functional [26]. The con¬ 
figurational entropy is estimated using the mean-field 
approximation [27]. Short-range clustering effects in 
Ti 1 _ a ,Al x N have been shown to be significant and similar 
in magnitude to the mean field contribution at temper¬ 
atures below 3500 K [19]. In principle, such clustering 
effects should be included in the free energy calculations. 
However, in Tii_ x Al x N they were found to predomi¬ 
nantly impact the AlN-rich compositions, having little 
effect on the solubility limits in the TiN-rich region, and 
only moderately decrease the maximum temperature of 
the miscibility gap. For these reasons, short-range clus¬ 
tering are not included in our calculations. 

We decompose the Helmholtz free energy into vibra¬ 
tional (“vib”), configurational (“config”) and electronic 
terms (“el”): 

F = U — TS = F e i + F v ib + Econfig (3) 

Using the force constants in Eq. 2, we can calculate the 
phonon density of states g(uj, T), from which we compute 
the vibrational contribution to the Helmholtz free energy: 


r OO fc 

F vih (u,T) = AU 0 (T) + g(u},T)^-du]+ 

( 4 ) 

+ / g{u,T)k B Tln ^1 - exp djjJ - 

In our case g(ui, T ) is strongly temperature dependent 
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and includes anharmonic effects and thermal expansion. 
AU 0 is an anharmonic correction to the energy in TDEP 
formalism that is added to (see supplementary ma¬ 
terial). 

The first test of our model system was to see that at low 
temperatures it does not deviate from a direct harmonic 
calculation using standard methods [9, 10], as illustrated 
in Fig. 1. The vibrational free energy of our model system 
agrees well with the harmonic model, with a difference of 
2 meV/atom at 0 K and 4 meV/atom at 1500 K. For this 
to be a fair test, we used the TDEP force constants from 
300 K where anharmonic contributions are expected to 
be negligible. This was a probe of the error introduced 
by the imposed symmetry, which turned out to be small. 

At higher temperatures, the harmonic model is no 
longer a good benchmark, and we calculate the formally 
exact free energy difference using thermodynamic inte¬ 
gration (UP-TILD) [24], Due to its considerable compu¬ 
tational cost, UP-TILD was applied to a single point, for 
which anharmonic effects are considerable. The correc¬ 
tion to the free energy (Eq. 1) at 1500 K was calculated 
to be ~ 5 meV/atom. We note that the harmonic ap¬ 
proximation is ~ 15 meV/atom further from the formally 
exact number (See Fig. 1 in supplementary material). 
We conclude that the SIFC-TDEP and UP-TILD free 
energies are within our desired error bars, and can use 
our model system free energies for the subsequent calcu¬ 
lations with confidence. 

We construct the Gibbs free energy (G = F+PV) from 
the Helmholtz free energy surface, where P = — dF/dV. 
The Gibbs free energy of mixing is given by, 

x) = Gxi 1 _ a .Ai x N — (1 — £)GTiN - zGain, (5) 

where x is the fraction of AIN. 

We reconstruct the concentration-temperature phase 
diagram at zero pressure by a direct minimization of G m ; x 
at each temperature and global composition on a grid 
(see supplementary material). We determine the spin- 
odal from the condition ( 9 |t < 0). 

The calculated phase diagram is compared with one 
neglecting lattice vibrations in Fig. 2. The difference 
in both the spinodal region and the miscibility gap is 
dramatic. By including the vibrational contribution to 
the free energy, the maximum of the miscibility gap is 
lowered from 6560 K to 2860 K. The roots of the second 
derivative of the total Gibbs free energy below 1000 K 
depend on the fitting function, in this case third order 
Redlich-Kister polynomials [28] (gray lines of spinodal). 
t 

Experimental verification of theoretical phase diagram 
(Fig. 2) is a challenging task. Because of the high melting 
temperature of TiN, its alloys are synthesized by thin- 
film deposition techniques, often followed by an anneal¬ 
ing. Due to the explicit out-of-equilibrium nature of this 
metastable transformation, direct information from equi¬ 
librium experiments is not available, so we adopted the 
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FIG. 2. (Color online) Calculated phase diagram for Bl- 
Tii_ x AUN. Comparing the dotted and solid lines shows the 
effect of including the vibrational entropy for binodal using 
the SIFC-TDEP method. The dash-dotted line corresponds 
to the spinodal metastable region including the anharmonic 
contribution, and the dashed line in solid solution region to 
the spinodal line without the vibrational contribution. 
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FIG. 3. (Color online) The limits of the concentration of 
metal sublattices as a function of aging time. 


following strategy. 4 /im thick TiAIN films with equal 
amounts of Ti and A1 (x = 0.5) were deposited by ca¬ 
thodic arc evaporation, resulting in a nearly uniform solid 
solution [29]. Our theoretical predictions suggest that 
when diffusion is activated this alloy will phase segregate 
through spinodal decomposition. To study this effect, 
a set of samples was prepared by annealing the TiAlN- 
films at 1073 K for 2, 4, 12 and 100 hours in an argon 
atmosphere. The evolved microstructures where then ex¬ 
amined by atom probe tomography. 

Fig. ?? shows 2D compositional maps of 2 nm thick 
slices through the reconstructed 3D atom probe speci¬ 
men of the annealed samples. The initially homogeneous 
cubic solid solution decomposes isostructurally forming 
cubic Ti- and Al-rich domains during the first 12 hours. 
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After 100 hours, the metastable cubic AlN-rich phase has 
to a large extent transformed to its stable wurtzite phase 

[21]. During the first 12 hours the evolution of the mi¬ 
crostructure and the chemical segregation becomes grad¬ 
ually more pronounced, which is consistent with the spin- 
odal decomposition [30]. After 100 hours the microstruc¬ 
ture is coarser. 

The phase diagram (Fig. 2) shows a substantial shift 
in the solubility at 1073 K compared to the predictions 
excluding vibrational effects, especially for the solubil¬ 
ity of A1 in TiN. In order to experimentally verify this, 
we extracted the maximum and minimum local compo¬ 
sition in the decomposed microstructure by constructing 
histograms of the concentrations in 1 nm 3 bins through¬ 
out the sample. These two local compositional extrema 
are plotted in Fig. 3 as a function of annealing time, 
and demonstrate the timescale and decomposition path 
of the alloy. The solubility of A1 in TiN at 1073 K tends 
to a final solubility of ~2 % according to the theoreti¬ 
cal phase diagram. The solubility limit for Ti in cubic 
AIN could not be determined experimentally since it was 
not reached prior to the onset of wurtzite AIN forma¬ 
tion. However, the experimental data for Ti in cubic AIN 
asymptotically approaches the theoretical value (x = 1). 

In summary, we present an accurate technique for cal¬ 
culating the vibrational contribution to the Gibbs free 
energy for random alloys. We reconstruct the phase dia¬ 
gram of a model Tq^ARN alloy, and demonstrate that 
in this system, the non-harmonic vibrational phonon free 
energy is large and comparable to the harmonic energies. 
As a result, we find a dramatic decrease of the maxi¬ 
mum temperature for the miscibility gap, from 6560 K 
to 2860 K, as well as an increase AIN of solubility in TiN 
as compared to calculations which neglect lattice vibra¬ 
tions. Atom probe tomography experiments on annealed 
Tio.5Alo.5N samples are in line with our theoretical pre¬ 
dictions demonstrating a finite AIN content in the TiN 
rich compositions after 100 h of annealing at 1073 K. 
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